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SCATTERING  BY  PERIODIC  SURFACES  WITH  SINUSOIDAL 
HEIGHT  PROFILE-A  THEORETICAL  APPROACH 

PART  II:  Numerical  Evaluation 

1.  INTRODUCTION 

In  part  I of  this  study,  a theoretical  approach  to  the  problem  of 
scattering  by  periodically  corrugated  metal  surfaces  is  presented  which  is 
basically  an  extension  of  the  physical  optics  approximation  leading  to  a 
rigorous  theory.  The  theory  assumes  the  current  distribution  in  the  metal 
surface  to  be  the  physical  optics  current  density  modified  by  multiplica- 
tion with  a Fourier  series  whose  fundamental  period  is  equal  to  the  period 
of  the  surfaces  grooves  (Floquet's  theorem).  The  coefficients  of  the 
Fourier  series  are  determined  from  the  extended  boundary  condition  which 
requires  that  the  field  radiated  by  the  current  distribution  into  the  half- 
space below  the  metal  surface  cancel  the  incident  field  in  this  half-space. 

The  method  reduces  the  scatter  problem  to  the  familiar  problem  of  solving 
a linear  system.  A particular  advantage  of  this  method  is  that  for 
certain  basic  types  of  height  profiles,  including  the  sinusodial  profile 
considered  here,  all  matrix  coefficients  of  the  linear  system  reduce  to  closed 
form  expressions  of  well-known  functions.  The  theory  is  thus  amenable  to 
efficient  computer  evaluation. 

A computer  program  has  been  written  [1]  and  is  operational.  Numerical 
evaluations  have  been  performed  on  the  ECOM  Burrough's  B-5500  computer. 

Program  outputs  include  the  current  distribution  in  the  metal  surface,  the 
complex  amplitudes  of  the  plane  wave  spectrum  (space  harmonics  representa- 
tion) of  the  scatter  field  and  an  accuracy  check  based  on  conservation  of 
power.  This  data  has  been  computed  for  both  TE-andTM-polarization,  where 
TE  and  TM  is  understood  with  respect  to  the  direction  of  the  surface  grooves. 

The  study  has  been  performed  in  support  of  work  on  a practical  problem, 
the  polarization  question  of  microwave  landing  systems.  One  of  the  aspects 
of  this  problem  concerns  scattering  by  large  periodic  metal  surfaces,  such 
as  hangar  doors,  and  the  question  whether  the  specular  reflection  coef- 
ficient of  such  surfaces  can  be  significantly  reduced  by  appropriate  choice 
of  polarization.  The  numerical  results  confirm  experimental  evidence 


The  program  is  based  on  a formulation  of  the  scatter  problem  using  a 
Fourier  representation  of  the  current  distribution  in  terms  of  sine 
and  cosine  functions  rather  than  in  terms  of  the  complex  exponential 
functions  used  in  the  theory  presented  in  Part  I.  The  latter  re- 
presentation leads  to  an  analytically  more  elegant  formulation  and, 
therefore,  has  been  preferred  in  presenting  the  theory. 


C2]  that  horizontal  polarization  generally  leads  to  substantially  less 
specular  reflection  than  vertical  polarization,  in  particular  in  the  interest- 
ing range  of  large  angles  of  incidence  (near  grazing).  It  is  implied  here  that 
the  grooves  of  the  periodic  surface  run  in  the  vertical  direction.  In  the 
limiting  case  of  grazing  Incidence,  the  specular  (amplitude)  reflection  co- 
efficient of  course  approaches  -1  for  both  polarizations. 

Additional  numerical  data  can  be  found  in  a preceedlng  summary  report  [3] 
which  also  discusses  the  case  of  circular  polarization  of  the  incident  field. 

It  is  shown  that  this  polarization  is  highly  effective  in  suppressing  higher 
order  grating  lobes  having  the  same  polarization  (circular  with  the  same  sense 
of  rotation)  as  the  incident  wave.  Grazing  lobes  with  the  opposite  sense  of 
rotation  are,  in  general,  present  with  significant  amplitudes.  Such  grazing 
lobes,  however,  would  not  be  received  by  an  airborne  antenna  polarized  for  opti- 
mum reception  of  the  primary  microwave  beam. 

Throughout  Part  II,  the  notation  of  Part  I is  used.  To  illustrate  the 
scatter  problem  considered  and  the  coordinates  used.  Fig.  1 of  Part  I is  re- 
produced. As  in  Part  I,  it  is  assumed  that  the  profile  of  the  metal  surface 
has  the  form 

z =*  z (x)  = h sin(27T  ^) 

a 

and  that  the  primary  field  is  a plane  wave  of  (suppressed)  time  dependance 

i(l)t  » «*-••  '*  * 

e incident  in  a direction  parallel  to  the  x,z-*plane.  The  scatter  problem 
for  other  directions  of  incidence  can  be  mathematically  reduced  to  this  case. 

2.  CURRENT  DISTRIBUTION 

In  Figures  2a  through  9b,  the  current  distribution  over  one  period  of  the 
surface  is  plotted  as  a function  of  the  lateral  coordinate  x.  The  curves 
apply  to  a surface  profile  with  parameters  d/X  ■ 2.5,  h/X  * 0.375  and  to  inci- 
dence angles  0=0°,  30°,  60°,  and  85°.  The  current  distribution  is  periodic 
with  the  period  d of  the  surface,  except  for  a phase  factor,  exp  (-ikx  sin  0), 
determined  by  the  incident  plane  wave  (Floquet's  theorem):  see  eqs.  (15a)  to 
(15c)  of  Part  I.  In  the  case  of  TE-polarization  where  the  incident  electric 
field  strength  has  been  assumed  to  have  unit  amplitude,  the  modulus  of  the 
current  density  is  normalized  by  multiplication  with/y0/e0;  in  the  case  of 
the  TM-polarization,  the  incident  magnetic  field  strength  has  been  assumed 
to  have  unit  amplitude  and  no  normalization  is  required. 


[2]  P.  S.  Demko,  "Polarization/Multipath  Study,  Technical  Memorandum 
VL-5-72,  Avionics  Laboratory,  U.  S.  Army  Electronics  Command,  Fort 
Monmouth,  New  Jersey,  1972. 

[3]  F.  Schwering  and  G.  Whitman,  "a  Theory  of  Scattering  by  Sinusoidal 
Metal  Surfaces",  Tech.  Rept.  ECOM-4496,  U.  S.  Army  Electronics 
Command,  Fort  Monmouth,  N.  J.,  May  1977. 
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In  the  valleys  of  the  surface  corrugations,  i.  e.,  in  the  range  - 4 < x < 0, 
TE-polarization  leads  to  appreciaoly  smaller  and  less  varying  current  2 
densities  than  TM— polarization.  This  effect  is  pronounced  in  particular  near 
normal  incidence.  Near  grazing  incidence,  the  surface  currents  for  both  TE  - 
and  TM  - polarization  are  substantially  reduced  in  magnitude  over  the  entire 
metal  surface. 

3.  SPACE  HARMONICS 

a.  Power  vs  Incidence  Angle 

In  Figures  11a  through  12b  and  14a  through  15b,  the  powers  P of  the  propa- 
gating space  harmonics  of  the  scatter  field  are  plotted  as  functions  of  the 
incidence  angle  8 of  the  primary  wave.  Two  surface  profiles  are  considered; 
their  parameters  are  d/X  **  2.5,  h/X  = 0.375  and  d/X  * 1.3,  h/X  = 0.1333,  re- 
spectively. P in  these  figures  denotes  the  power  transmitted  by  the  mth  space 
harmonic  through  unit  area  of  planes  z =»  const,  assuming  that  the  incident 
power  per  unit  area  is  unity.  Hence  in  the  notation  of  Part  I (see  eqs.  (20a) 
and  (20b)): 


t 

E (D|2 

m j cos6m 

for  TE-polarization 

(la) 

p = < 

m 1 

cos8 

< 

H(1)l2 

m | cos8m 

for  TM-polarization 

(lb) 

cos  8 


The  direction  angles  8m  of  the  propagating  space  harmonics  versus  incidence 
angle  8 are  shown  in  Figs.  10  and  13  for  the  two  surface  profiles  considered. 
These  curves  apply  to  both  TE-and  TM-polarization. 

Figures  11a  to  12b  and  14a  to  15b  show  that  in  the  range  of  large  incidence 
angles  (8  > 60°)  the  specular  reflection  coefficient  (P0)  is  significantly 
greater  for  TE-than  for  TM-polarization.  This  result  is  significant  in  the  con- 
text of  the  polarization  problem  of  microwave  landing  systems  mentioned  in  the 
Introduction.  Thus  undesired  specular  reflection  of  microwave  beams  from 
periodic  structures,  such  as  hanger  doors  at  airports,  can  be  substantially  re- 
duced by  appropriate  choice  of  polarization.  The  surface  parameters  used  in 
the  numerical  evaluation  approximate  those  of  actual  profiles* 

b . Anomalies 

The  power  distribution  curves  Pm  vs  0,  show  Rayleigh-type  Wood  anomalies 
as  discussed,  for  example,  by  Hessel  and  Oliner  [4].  At  angles  of  incidence, 
where  one  of  the  space  harmonics  is  scattered  Into  the  direction  of  grazing 
and  thus  changes  from  being  a propagating  wave  to  being  an  evanescent  one,  or 
vice  versa,  the  power  distribution  over  the  remaining  space  harmonics  is  radi- 


[4]  A.  Hessel  and  A.  A.  Oliner,  "A  New  Theory  of  Wood  Anomalies  on  Optical 
Gratings"  Applied  Optics,  4,  pp.  1275-1297,  1965. 
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cally  altered  and  one  or  several  spectral  order  may  attain  very  high  or 
near  zero  Intensities.  The  incidence  angles  at  which  such  anomalies  occur 
are  found  from  the  relation  for  the  propagation  angles  of  the  space  har- 
monics (eq.  (3a),  Part  I) 

sin  0 * sin  0 + m -3- 

m d 

by  letting  ©m  = +90°.  Thus 

0 “ 0^  '■  -arc  sin  (m  — + 1) 

These  angles  are  indicated  in  Figs.  10  through  15b  by  small  arrows.  An 
additional  anomaly  occurs  at  0 = 90°  where  the  specular  reflected  wave  is 
at  grazing  angles. 

The  anomalies  are  pronounced  for  TM-polarization  (S-type  anomalies)  and 
weaker  for  TE-polarization  (P-type  anomalies ).  This  can  be  explained  by 
observing  that  - in  particular,  w]ien  b«X  - the  electric  field  of  a grazing 
TE  wave  is  practically  short  circuited  by  the  nearly  plane  conducting  sur- 
face [5] . These  results  are  in  agreement  with  the  experimentally  well 
established  fact  that  P-type  anomalies  require  substantially  deeper  surface 
corrugation  than  S-type  anomalies  to  become  observable.  For  large  d/X-  values, 
Rayleigh  type  anomalies  tend  to  become  less  pronounced,  which  is  probably  due 
to  the  larger  number  of  space  harmonics  involved  in  the  redistribution  of 
power.  Note  that  for  d/X  =2.5  both  anomalies  are  double  anomalies  [6]  (two 
different  space  harmonics  are  at  grazing  angles  simultaneously)  so  that-  a — 
particularly  strong  effect  can  be  expected.  For  d/X  = 1.3  only  single  anom- 
alies occur . 

A second  type  anomally  expected  to  occur  is  of  the  "Brewster  angle"  type 
[7] ; it  can  be  characterized  as  a resonance  between  the  incident  wave  and  one 
of  the  space  harmonics  [51.  The  condition  for  this  type  of  resonance  to  take 
place  is  that 

0 = 0„  = arc  sin  {(n+%)  -j  } (3) 

D U , 


( 51  A.  Hessel ^ and  J.  Shmoys,  "Computer  Analysis  of  Propagation  Reflection 
Phenomena ^Polytechnic  Institute  of  Brooklyn  Scientific  Report  Contract 
Number  AAB07-73-19-2716,  August  1973. 


[6]  Double  anomalies  occur  when  d/X  is  an  integer  or  half  integer;  see  eq. 


[7]  D.  Y.  Tseng,  A.  Hessel,  and  A.  A.  Oliner,  "Scattering  by  a Multimode 
Corrugated  Structure  with  Applications  to  P-type  Wood  Anomalies?  Alta 
Frequenza,  38,  1968  Special  Issue  on  the  URSI  Symposium  on  Electro- 
magnetic Waves,  Stress,  Italy. 


(2a) 

(2b) 
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but  it  also  depends  on  the  depth  os  the  surface  corrugations.  It  may  occur 
for  both  TE-  and  TM-polarization  but  has  been  found  to  be  more  sensitively 
dependent  on  h in  the  latter  case  [8]. 

For  d/X  » 2-1’,  conditions  (2)  and  (3)  yield  the  same  values  for  0^  and 
Qg  so  that  Rayleign  and  Brewster  angle  anomalies  may  reinforce  each  other  but 
are  difficult  to  distinguish.  For  d/X  « 1.3  the  Brewster  angle  condition 
yields 

0B  - 22.52° 

Observe  the  decline  of  P (the  specular  reflected  power)  and  the  increase  of 
P_1  near  this  angle  in  the  case  of  TM-polarization;  see  Figs.  15a  and  15b. 
Furthermore,  Figs.  19a  and  19b  show  that  if  the  surface  depth  is  appropriately 
chosen,  i.e.  h/X  ■ 0.19  for  TM-  and  h/X  = 0.27  for  TE-  polarization,  the 
specular  reflected  power  becomes  zero  at  this  angle  so  that  all  power  is  pro- 
pagated by  the  space  harmonic  of  order  m = -1.  This  means  that  the  Incident 
power  is  backscattered  into  the  direction  of  arrival.  From  eqs.  (2a)  and  (3) 
it  can  be  seen  that  0_^  ■ -0  (backscatter  condition)  when  the  angle  of  arrival 
coinsides  with  the  zero  order  Brewster  angle,  the  situation  considered  in 
Figs.  19a  and  19b. 

c.  Dependence  on  Surface  Depth 

Figs.  16a  to  19b  show  the  power  propagated  by  the  various  space  harmonics 
of  the  scatter  field  as  functions  of  the  depth  of  the  surface  grooves.  Para- 
meters are:  d/X  = 2.5;  0 * 30°,  60°,  85°  and  d/X  ■*  1.3;  0 « 22.62°  (Brewster 
angle).  Note  that  the  groove  depth  is  2h  rather  than  h,  the  parameter  used  in 
the  figures. 

For  both  polarizations,  the  dependence  of  P on  h exhibits  a basically 
oscillatory  behavior,  though  TM-polarization  tencfs  to  be  more  sensitive,  in 
particular  at  angles  of  incidence  near  grazing;  see  Figs.  18a  and  18b.  These 
effects  have  been  discussed  by  Zaki  and  Neureuther  [91. 

When  the  incidence  wave  is  arriving  at  the  Brewster  angle,  (Figs.  19a  and 
19b)  the  power  of  the  specular  reflected  wave,  PQ(h),  has  zeros  for  both 
polarizations,  as  has  been  discussed  above.  However,  in  the  TE-  case,  deeper 


[8]  See  reference  [5]  on  page  4. 

[9]  K.  A.  Zaki  and  A.  R.  Neureuther,  "Scattering  From  a Perfectly  Conduct- 
ing Surface  with  Sinusodial  Height  Profile, TM-Polarization>IEEE  Trans. 
Antennas  and  Prop.,  AP-19,  pp.  747-751. 
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grooves  ere  required  to  ettsln  a zero  than  in  the  TM  - case.  Similarly,  Figs. 
16a  and  16b  show  near-zeros  occurring  at  h/X  “ 0.19  for  TM-  and  h/X  - 0.26 
for  TE-polarizatlonUO]  . These  near-  zeros  apparently  are  due  to  the  inci- 
dence angle,  6 “ 30°,  being  rather  close  to  the  Brewster  angle  at  8 - 36.87°. 
Note  that  TM-polarization  leads  to  a near  - zero  also  at  8 - 60°  though  no 
Rayleigh- type  or  Brewster  angle  anomalies  are  expected  in  the  viciaty  of  this 
angle.  This  near  - zero  may  be  attributed  to  a broad  minimum  of  PQ  (as  a func- 
tion of  8)  In  the  range  of  large  incidence  angles;  see  Fig.  12a.  v 

4.  ACCURACY  CHECK 

Two  criteria  are  used  to  check  on  the  accuracy  of  computed  data:  conserva- 
tion of  power  and  reciprocity. 

a.  Conservation  of  Power 

Since  the  metal  surface  is  assumed  perfectly  conducting,  the  total  scatter 
power,  i.e. , the  accumulated  power  of  all  propagating  space  harmonics  must  equal 
the  power  of  the  incident  wave.  Hence,  with  eqs.  (la)  and  (lb) 

!CPm  5 cose  XjEtt<X>  ^ cosem  " 1 (4a) 

m m 

for  TE-polarization  and 

$ ^ Dh.<1)|2  1 <«>> 

m m 

for  TM-polarization.  The  summation  is  taken  over  the  propagating  spectral 
orders,  i.e.,  over  all  Integer  m within  the  range 

- ^ (1  + sin  8)  < m < ^ (1  - sin  8). 

See  eq.  (3b)  of  Part  I. 

In  Tables  1 and  2,  the  relative  power  errors 

^ Y lE  (1>I2  cos6m  - 1 (3a) 

TE  cos8  4_«  1 m 1 m 

m 

and 

eiM  ■ £ |«,(1)|2  “»e„  - 1 (3b> 

® m 


no]  It  is  interesting  to  note  that  the  groove  depths  h at  which  these  zeros 
or  near  zeros  occur  are  practically  the  same  for  d/X  ■ 1.3  and  2.5,  the 
two  cases  considered  here. 
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are  listed  as  functions  of  6 for  (1)  d/X  = 2.5,  h/X  = 0.375  and  (2) 
d/X  * 1.3,  h/X  = 0.1333,  i.e.,  for  the  surface  parameters  considered  in 
most  of  the  figures.  In  addition  to  the  power  errors,  the  highest  peaks 
of  the  current  distributions  in  the  metal  surface  are  shown.  Near  the 
angles  6 * 9A,  where  Rayleigh- type  Wood  anomalies  are  expected  to  occur, 
smaller  increments  in  8 are  chosen  in  the  tables  to  make  effects  apparent. 

As  indicated  in  the  figures,  we  have  two  Rayleigh  wavelength  angles  for 
d/X  * 2.5  i.e. , 

0A  - 11.54°  and  36.87° 

However,  since  d/X  is  a half  integer,  the  associated  anomalies  are  de- 
generate (double  anomalies).  For  d/X  - 1.3  only  single  anomalies  occur; 
the  Rayleigh  wavelength  angles  are 

8A  - 13.34°  and  32.58° 

The  tables  show  that  the  accuracy  attained  with  the  new  method  [11]  is 
excellent  for  parameters  (2),  and  acceptable  for  parameters  (1)  apart  from 
a small  range  of  incidence  angles  near  the  second  anomaly  where  the  power 
error  in  the  TM  case  reaches  values  in  the  order  of  10%.  In  general,  the 
accuracy  for  TE-polarization  is  higher  than  for  TM-polarization,  at  most 
incidence  angles  by  an  order  of  magnitude.  Note  that  passing  an  anomaly  is 
usually  accompanied  by  a sharp  drop  or  increase  in  accuracy  while  the  current 
amplitude  tends  to  peak  up  near  anomalies;  see  in  particular.  Table  1.  The 
drastically  increased  current  amplitude  shown  in  this  Table,  at  the  second 
anomaly  however,  may  not  be  real.  In  several  instances  it  has  been  observed 
that  a beginning  loss  of  accuracy  expresses  itself  more  sensitively  in  a 
rapidly  increasing  amplitude  of  the  current  distribution  rather  than  in  an 
increasing  violation  of  the  power  criterion,  an*  effect  which  commences  only 
later.  As  pointed  out  already  in  the  discussion  of  Figs.  8a  and  9a,  the 
current  amplitudes  decline  substantially  for  both  polarizations  as  0 ap- 
proaches 90°. 

b.  Reciprocity  Theorem: 

Consider  two  plane  waves  incident  from  directions  0a  and  0b,  respectively. 
If  the  first  of  these  plane  waves  scatters  a space  harmonic  of  order  m into 
the  direction  of  the  second  wave  such  that  [12]  0®  - - 0b,  then  the  second 


[11]  The  power  errors Cjg  and  depend  , of  course,  on  the  degree  to  which 
the  linear  systems  (18c)  and  (18d)  of  Part  1 are  truncated  in  the 
numerical  evaluation.  The  values  listed  in  Tables  1 and  2 refer  to 
systems  of  sizes  25  x 25  and  22  x 22,  respectively. 

[12]  For  explanation  of  the  minus  signs  in  the  relations  0a  — 0b  and  6b  =-8a 
we  refer  to  Fig.  1.  Incidence  angles  of  the  primary  field  are  counted 
positive  in  the  2nd  quadrant  of  the  x,z-plane  while  the  direction  angles 
of  the  scatter  field  are  counted  positive  in  the  1st  quadrant. 
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plane  wave  produces  a space  harmonic  of  the  same  prder  m propagating  in 
the  direction  of  the  first  incident  wave,  hence  0”  - -0a,  and  the  complex 
amplitudes  of  the  two  space  harmonics  are  related  by: 


(Aa) 


for  TE-polarization  and 


(Ab) 


for  TM-polarizatlon.  It  is  assumed  here  that  the  two  incident  waves  have 
equal  amplitude  and  phase.  A specific  advantage  of  the  reciprocity  rela- 
tions is  that  they  check  amplitude  and  phase  of  individual  space  harmonics. 

In  Table  3 and  A the  results  of  a reciprocity  check  are  shown;  surface 
parameters  are  d/X  - 2.5,  h/X  - 0.375.  In  the  TE-  case,  relative  amplitude 
errors  are  generally  below  0.5%  - larger  errors  are  due  to  small  amplitudes- 
and  phase  errors  remain  below  1°.  In  the  TM-  case  discrepancies  are  ap- 
preciably larger  in  conformance  with  a trend  observed  when  applying  the  power 
criterion.  In  general,  it  has  been  found  that  when  the  power  criteria  is 
well  satisfied  the  same  is  true  for  the  reciprocity  criterion,  though  specific 
power  and  reciprocity  errors  do  not  (and  are  not  expected  to)  directly 
correlate. 


c.  Limitations  of  Method: 

The  accuracy  of  the  computed  results  has  been  found  to  be  sensitively 
dependent  on  the  depth  of  the  surface  grooves  i.e.,  on  the  parameter  h/X. 

The  power  and  reciprocity  criteria  are  excellently  satisfied  at  small 
h/X  <0.1  and  well  satisfied  up  to  h/X  ^0.5.  However  when  the  surface 
depth  approaches  and  is  Increased  beyond  one  wavelength,  the  accuracy  sooner 
or  later  breaks  down  (depending  on  d/X  and  9).  This  failure  occurs  for  TE- 
and  TM-  polarization  at  approximately  the  same  h/X  values.  To  demonstrate  the 
decline  of  accuracy  with  increasing  h,  the  power  errors  •nd^TH  encountered 
in  computing  the  curves  Pm  - Pa(h)  of  Figs.  16a,  16b  and  I9a,  19b  are  listed 
in  Tables  5 and  6 as  functions  of  h/X.  It  has  been  found  that  the  present 
program  does  not  permit  solution  of  the  problem  of  declining  accuracy  by  the 
obvious  method  of  arbitrarily  increasing  the  si*e  of  the  linear  systems,  eqs. 
(18c)  and  (18d)  of  Part  I,  in  the  numerical  evaluation.  Up  to  now  it  has  not 
been  established  whether  the  accuracy  problem  is  due  to  Inherent  poor  con- 
vergence of  the  linear  systems  at  large  h/X  or  to  a break  down  of  aubroutinea 
employed  In  the  computer  program  such  as  those  for  computing  Bessel  functions 
of  complex  arguments.  Further  work  is  required  to  answer  this  question  and, 
if  possible,  to  extend  the  h/X  range  which  can  be  handled  by  the  program. 

No  limitations  have  been  encountered  with  regard  to  admissible  d/X  valuas. 
The  accuracy  as  measured  by  the  power  criterion  tends  to  increase  with  in- 
creasing d/X  which  would  be  expected  since  the  underlying  method  is  basically 
an  extension  of  the  physical  optics  approximation. 
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5. 


COMPARISON  WITH  EXISTING  METHODS 


Tables  7a  through  8d  show  a comparison  of  numerical  values  obtained  by  use 
of  the  new  method  (WS) , the  method  by  Tong  and  Senior  (TS),  and  the  physical 
optics  approximation  (PO).  These  tables  have  been  taken  directly  from  Tong 
and  Senior  [13 ] where  only  the  values  computed  according  to  the  new  technique 
have  been  added  [14 ] . Tables  9 and  10  compare  the  accuracy  of  these  approaches 
in  satisfying  the  power  criterion. 

The  tables  show  that  the  WS  and  TS  methods  in  general  are  in  good  agree- 
ment (apart  form  apparent  sign  errors  in  a number  of  cases),  while  the  PO 
approximation  leads  to  substantial  errors  as  has  been  pointed  out  already  by 
Tong  and  Senior.  Such  errors,  however,  are  not  unexpected  since  the  periods 
d of  the  surfaces  considered  here  do  not  satisfy  the  condition  d » X under 
which  application  of  the  PO  approxlmtion  would  be  justified.  Taking  shadow 
effects  into  account  appears  to  increase,  rather  than  decrease  errors  at  the 
small  d/X  ratios  considered  here.  Tables  9 and  10  show  that  the  accuracy  of 
the  WS  method  in  general  exceeds  that  of  the  TS  method  though  improvements 
are  not  dramatic.  There  are  exceptions;  for  example,  for  d/X  - 0.4,  h/X  - 0.2, 
0 ■ 0°,  the  TS  method  satisfied  the  power  criterion  with  smaller  error  than 
the  WS  method.  On  the  other  hand,  for  d/X  - 0.2,  h/X  ■ 0.1,  0 ■ 60°,  the  WS 
method  exceeds  the  accuracy  obtained  with  the  TS  method  by  an  order  of  magni- 
tude. For  a comparison  of  the  WS  program  with  the  program  developed  by  7,aki 
and  Neureuther  see  reference  [3]. 

Tables  11a  through  lid  test  the  accuracy  of  the  PO  approximation  in  the 
case  of  a surface  satisfying  the  condition  d » X.  Parameters  are 

| - 10;  | - 0.5;  0 - 15°,  45°,  75°,  and  85°. 

Listed  are  the  power  densities  Pm  of  the  propagating  space  harmonics  com- 
puted according  to  the  WS  and  P0[15]  methods,  respectively.  For  0 - 75°  and 


[13]  T.  C-H.  Tong  and  T.  B.  A.  Senior,  "Scattering  Electromagnetic  Waves  by 
a Periodic  Surface  with  Arbitrary  Profile,  Scientific  Report  AFCRL-72- 
0258,  April  1972,  prepared  by  the  University  of  Michigan,  Radiation 
Laboratory. 

[141  Tong  and  Senior  assume  a cos-profile  of  the  metal  surface;  the  present 
theory  uses  a sin-profile.  Transition  form  a sin-to  a cos-profile  is 
accomplished  by  multiplication  of  the  complex  amplitudes  of  the  space 
harmonics  by  a phase  factor  exp  (-lmir/2).  This  adjustment  has  been  made 
in  Tables  7a  through  8d  for  the  WS  values.  In  the  TE-  case  however,  a 
discrepancy  by  180°  remains  which  can  probably  be  traced  to  a sign  error. 

[15]  For  a summary  of  the  P0  approximation  method  see  Appendix  B of  Part  I. 
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85°,  geometrical  shadowing  occurs  and  PO  values  are  given  for  both  options, 
that  shadowing  is  considered  and  that  it  is  neglected.  Only  power  densities 
exceeding  10  3 (referred  to  the  incident  power  density)  are  shown  in  the 
tables;  zeros  indicate  smaller  magnitudes.  Due  to  the  large  d/X  ratio  the 
number  of  propagating  space  harmonics  is  rather  large  (20).  At  the  bottom 
of  the  tables,  the  accumulated  power  of  all  space  harmonics  and  the  magni- 
tude of  the  relative  power  error 


iei-|Ep„-  x|  • 1001 

are  given. 

The  tables  show  that  for  directions  of  incidence  sufficiently  far  above 
grazing  the  PO  method  provides  a good  approximation  of  the  power  densities 
of  the  space  harmonics,  though  the  accuracy  declines  with  increasing  0.  At 
0 * 45°  and  75°,  errors  of  10%  or  more  occur  even  for  space  harmonics  of 
significant  magnitude.  Near  grazing  incidence,  at  0 • 85°,  the  PO  approxima- 
tion with  shadowing  leads  to  substantial  errors,  but  may  still  be  used  to 
estimate  the  order  of  magnitude  of  the  space  harmonics.  The  PO  approximation 
without  shadowing  is  practically  useless  in  this  range;  the  power  criterion 
is  violated  with  an  error  exceeding  100%.  The  apparent  reason  for  this  failure 
is  that,  at  incidence  angles  near  grazing,  coupling  between  surface  cells 
becomes  significant,  in  particular  at  TM-polarizatlon  where  even  a slight 
surface  corrugation  has  the  effect  of  Inverting  the  (amplitude)  reflection 
coefficient  of  a perfectly  plane  metal  surface  from  +1  into  -1.  In  con- 
trast, the  PO  approach  is  based  on  the  assumption  that  the  current  distribu- 
tion is  determined  solely  by  local  effects.  In  general,  it  appears  that  tak- 
ing shadow  effects  into  account  improves  the  approximation  at  large  d/X. 

Note  that  for  all  incidence  angles  considered  here  the  deviation  of  the  total 
scatter  power  from  unity  (power  criterion)  is  smaller  than  the  power  errors 
of  Individual  space  harmonics.  Apparently  these  errors  compensate  to  a certain 
degree  in  the  accumulation  process. 

It  is  assumed  here  that  the  Pm  values  computed  on  the  basis  of  the  WS 
theory  are  sufficiently  accurate  to  serve  as  a standard  for  judging  the 
accuracy  of  the  PO  approximation.  As  justificaiton  we  point  to  the  very  high 
accuracy  with  which  the  WS  values  satisfy  the  power  criterion  at  all  0 values 
considered.  For  a further  discussion  of  the  limitations  of  the  PO  approxima- 
tion we  refer  to  Zakl  and  Neureuther  [6] 

Observe  that  at  large  angles  of  incidence  the  magnitude  of  the  specular 
reflected  wave  (spectral  order  m ■ 0)  is  very  small  for  TM-polarization,  but 
dominant  for  TE-polarization  [16] . This  is  in  agreement  with  the  trend, 
noted  for  surfaces  of  moderate  groove  width,  that  in  the  range  of  large  inci- 
dence angles  TM-polarization  leads  to  substantially  less  specular  reflection 
and  a more  even  distribution  of  the  Incident  power  over  the  spectrum  of  space 
harmonics  than  TE-polarization. 


[16]  For  0 » 90°,  the  power  of  the  specular  reflected  wave,  of  course, 
approaches  unity  for  either  polarization. 
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Accuracy  check  using  power  criterion.  Surface  parameters: 
d/X  » 2.5 1 “ 0.375.  In  addition  to  power  errors^  and^, 

table  shows  peak  amplitudes  of  current  distribution.  Ne«rRayleigh 
wavelength  angles  (0^  * ll.5^°  and  36. 87°), increments  In  0 are 
reduced  to  0.5°  to  make  effects  of  anomalies  apparent. 
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TABLE  2 

Accuracy  check  using  power  criterion.  Surface  parameters  d/\  - 1.3, 
h/\  * O.I333.  In  addition  to  power  errors  and  table  shows 
peak  amplitudes  of  current  distribution.  Near  Rayleigh  wavelength 
angles  (©A  * 13*34°  and  32. 58°) , increments  in  0 are  reduced  to  0.5°  to 
make  effectsof  anomalies  apparent. 
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Dependence  of  Power  Errors  Cmg  and  #5™  on  Surface  Depth 
Surface  Period:  d/ft  * 1*3 
Incidence  Angle:  t - 22.62°  CBrcVsfer 
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Comparison  of  approaches:  Oomplex  amplitudes  of  space 
harmonics.  (TE  - polarization;  d/3^0.2,  h/7*0.l) 
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Physical  Optics 
w/o  Shadow 
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0 
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Physical  Optics 
with!  Shadow 
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Comparison  of  approaches:  Oomplex  amplitudes  of  space 
harmonics.  (TE  - polarization;  d/Jl=l  , h/Xa0.25;  0=0°) 
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Comparison  of  approaches:  Oomplex  amplitudes  of  space  harmonics 
(TE  - polarization;  d/^0.4,  h/^r0.2) 
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TABLE  7d 

Comparison  of  approaches:  Oomplax  amplitudes  of  spaca  harmonics 
<TE  - polarization;  d/Jl*0.2,  h/a-0.03;  Mfi) 
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Ofemparison  of  approaches:  Oomplex  amplitudes  of  space 
harmonics.  (TM  - polarization;  d/^0.2,  h/JM).l) 


Mag.  Phase  Mag.  Phase 
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Comparison  of  approaches:  Oomplex  amplitudes  of  space 
harmonics.  (TM  - polarization;  d/^*l  f\,  h/*0.25;  0*0°) 
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I.  Whitman,  Schwering  0.9008  99.391  0.3320  +51 .93  1 0.3347 1 +-52.52 
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Cooper Ison  of  approaches:  Coop lax  amplitudes  of  space  harmonica  (TM  - 
polarization;  d/A  s 0*4,  h/‘A  ■ 0.2) 
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TABLE  8d 

Comparison  of  approaches:  Complex  amplitudes  of  space  harmonics  (TM  - 
polarization:  d/A  “ 0.2,  h/A  ■ 0.03;  • » 0°) 
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Scattering  from  surfaces  of  large  surface  period:  Compari 
of  W~S  theory  and  physical  optics  approximation. 

Surface  parameters:  d/^»IO,  h/^0.5 
Angle  of  incidence:  9= 1 5°. 


Spectral  Order 
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Power  Pm  of  Space  Harmonics 
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TABLE  lib 

Scattering  from  surfaces  of  large  surface  period:  Oomparison 
of  W-S  theory  artd  physical  optics  approximation. 

Surface  parameters:  d/^=IO,  h/^=0.5 
Angle  of  Incidence:  0=45° 


Spectral  Order 


Power  Pm  of  Space  Harmonics 


5P": 


W“S  Theory 


= 0.00  0.00 
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0.182  0.182 
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0.006  0.005 
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0.018  0.061 
0.333  0.225 


PO  Approximation 
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table  lid 

Scattering  from  surfaces  of  large  surface  period: 
Oomparison  of  W-S  theory  and  physical  optics  approximation. 
Surface  parameters:  d/Jl=IO,  h/3l=0.5 
Angle  of  Incidence: 


Spectral  Order 
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Power  Pm  of  Space  Harmonics 


pm=0.00 
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0.077 
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